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Mass transfer from a giant star to a main sequence 
companion and its contribution to long-orbital-period blue 
stragglers 
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ABSTRACT 

Binary population synthesis shows that mass transfer from a giant star to a main- 
sequence (MS) companion may account for some observed long-orbital period blue 
stragglers. However, little attention is paid to this blue straggler formation sce- 
nario as dynamical instability often happens when the mass donor is a giant star. In 
this paper, we have studied the critical mass ratio, q c , for dynamically stable mass 
transfer from a giant star to a MS companion using detailed evolution calculations. 
The results show that a more evolved star is generally less stable for Roche lobe over- 
flow. Meanwhile, q c almost linearly increases with the amount of the mass and angular 
momentum lost during mass transfer, but has little dependance on stellar wind. To 
conveniently use the result, we give a fit of q c as a function of the stellar radius at the 
onset of Roche lobe overflow and of the mass transfer efficiency during the Roche lobe 
overflow. 

To examine the formation of blue stragglers from mass transfer between giants 
and MS stars, we have performed Monte Carlo simulations with various q c . The 
simulations show that some binaries with the mass donor on the first giant 
branch may contribute to blue stragglers with q c obtained in this paper 
but will not from previous q c . Meanwhile, from our q , blue stragglers from 
the mass transfer between an AGB star and a MS companion may be more 
numerous and have a wider range of orbital periods than those from the 
other q c . 
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1 INTRODUCTION 

Blue stragglers (BSs) are stars that have remained on the 
main sequence for a time exceeding that expected, for their 
masses, from standard stellar evolution theory. The exis- 
tence of BSs indicates an incomplete understanding of stel- 
lar evolution and also of star formation within clusters 
(Stryker 1993). Since they are bright and blue, these ob- 
jects may affect the integrated spectra of their host clus- 
ters by contributing excess spectral energy in the blue 
and ultraviolet. BSs are therefore important in studies of 
population synthesis ( |Xin fc Deng 2005). Much evidence 
shows that these strange objects are relevant to primordial 
binaries l|Ferraro et al. 20031 |Davies, Piotto fe Angeli 2004| 
|MapelIi et al.2004| |. Binary coalescence from a contact bi- 
nary is a popular hypothesis for single BSs and it is be- 
lieved that these contact binaries are mainly from case A 
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mass transfer Q l|Mateo et al. 19901 IPols fc Marinus 19941 
lAndronov et al. 20061 IChen fc Han 2008[> . Mass transfer is 
another channel to produce BSs from primordial bina- 
ries. During Roche lobe overflow (RLOF), the less mas- 
sive star accretes some material and goes upward along the 
main sequence, if it is still a main-sequence star. The ac- 
creting component may then be observed as a BS when 
it is more massive than the turnoff of the host cluster 
(|Mc Crea 19641 IChen fc Han 20"04)) . Previous studies show 
that both case A and case B mass transfer are only respon- 
sible for some short- and mid-period BSs (Leona rd 19961 



1 According to the evolutionary state of the primary at 
the onset of mass transfer, three mass transfer cases 
are defined, i.e. case A for the primary being on the 
main sequence, case B for the primary after the main 
sequence but before central He burning, and case C 
for the primary during or after central He burning 
( Kippenh ahn &: Weigert 1967| ) . 
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Chen & Han 2004), which are rare in some old open clus- 
ters, and case C mass transfer may account for BSs in long- 
orbital period spectroscopic binaries. During RLOF, the or- 
bital period will decrease at first, however, if the primary 
continues to transfer material to the companion after the 
primary is less massive than the secondary, the binary will 
become wider and the orbital period will increase. Since dy- 
namical instability often happens when the mass donor is a 
giant star, there is little work to focus on for this channel to 
BSs. 

As is well known, a fully convective star will increase 
in radius with mass loss and decrease in Roche lobe ra- 
dius if it is more massive than its accreting companion 
(Paczynski 1965). This means that the mass donor will 
overfill its Roche lobe by an ever-increasing amount, leading 
to mass transfer on a dynamical time-scale, the formation 
of a common envelope (CE) and a spiral-in phase. The crit- 
ical mass ratio (the mass donor /the accretor) is about 2/3 
from a polytropic model with a polytropic index n = 1.5 
for conservative mass transfer, indicating that mass trans- 
fer would be dynamically unstable if the mass donor has a 
mass larger than 2/3 of the mass of the companion star. If 
stellar wind is not included before RLOF to decrease the 
primary's mass, case C mass transfer is always dynami- 
cally unstable for a binary with initial mass ratio qi > 1 
[3 (gi = M11/M21, Mu and M^ are the initial masses of the 
primary and of the secondary, respectively) , and therefore 
makes no contribution to BSs. However, there are some 
problems when we use the above criterion in a real binary 
system for the following facts. (1) Giant stars have large 
condensed cores (usually they are degenerate), so they can- 
not be modelled by a fully convective star. The critical mass 
ratio increases substantially because of the condensed core 
(Hjellming fc Webbink 1987] ). They conform to the formula 
as follows ( | Webbink 1988fl : 



q c = 0.362+ 



1 



3(1-Af /Afi) ' 



(1) 



where M c and Mi are the core mass and the total 
mass, respectively, of the donor as RLOF begins. (2) 

The condition for dynamical instability also depends on the 
amount of mass and angular momentum loss during RLOF. 
Assuming that the lost mass carries away the same specific 
angular momentum as pertains to the mass donor, Sober- 
man et al. (1997]) gave a fitted Roche lobe mass-radius expo- 
nent £l = —1.7/3 + (2.4/9 — 0.25)q, where /3 is mass transfer 
efficiency determined as the mass fraction of the lost mass 
from the primary accreted by the secondary. We may then 
obtain the critical mass ratio q c by setting the adiabatic 
mass-radius exponent £s = Cl> where £s can be fitted from 
the data of numerical calculations of Hjellming & Web- 
bink (1987) and Soberman et al. (|1997p (see also Han et al. 
2001). Fig.[l]clearly shows the dependence of q c on the mass 
transfer efficiency /3. (3) Stellar wind prior to RLOF, which 
will be strongly enhanced due to the tidal interaction with 
the companion, will increase the core fraction in equation 
(1) and allow the system to stabilize more easily. 

Possibly, a more fundamental problem of such a crite- 



2 The primary should be more massive than the secondary ini- 
tially from standard stellar evolution if we expect the binary to 
be composed of a giant star and a main-sequence companion. 




Figure 1. Critical mass ratio q c for dynamically unstable RLOF. 
/ c is the core mass fraction of the mass donor. The cross, plus, 
asterisk and circle lines are for mass transfer efficiency j3 = 0.25, 
0.5, 0.75 and 1.0, respectively. The lost mass is assumed to carry 
away the same specific angular momentum as pertains to the mass 
donor. The solid line is from Webbink (1988) (see equation (1), 
f c = M c /M±) for conservative RLOF. 



rion is that it does not take into account the detailed dy- 
namics of the mass transfer process. There are some sys- 
tems, observationally, which should experience dynam- 
ical mass transfer while appearing to avoid a CE phase 
( Podsiadlowski ct al. 1992). Several recent full binary evo- 
lution calculations also show that the simplistic criterion 
above is not really appropriate. For example, it is shown by 
Podsiadlowski et al. (2002) that mass transfer is dynamically 
stable for all giants up to a mass of about 2Mq , in the case 
of giants transferring mass to a neutron star of 1.3/1.4Mq. 
The study of Han et al. Q002) showed that q c may be up to 
about 1.3 as /3 = 0. 

In this paper, we will study the critical mass ratio for 
dynamical stability in mass transfer from a giant star to a 
main sequence companion from detailed binary evolutions, 
and adopt the results to estimate the contribution to BSs 
from the mass transfer. In section 2, we simply describe the 
numerical codes we have employed in the study. The re- 
sults of the critical mass ratio are shown in section 3. In 
section 4, we show some examples of binaries which even- 
tually form long orbital period BSs. The consequences from 
Monte Carlo simulations are shown in section 5. Discussions 
and conclusions are given in section 6. 



2 BINARY EVOLUTION CALCULATIONS 

The binary evolution code we employed was originally de- 
veloped by Eggleton ipW7T\ [T9T21 |T973) l. which has been 
updated with the latest input physics over the last three 
decades as described by Han et al. (|1994[> and Pols (|1995l 
I1998p . Some characteristics of the code, i.e. a self-adaptive 
non-Lagrangian mesh, the treatment of both the convec- 
tive and semi-convective mixing as diffusion processes, the 
simultaneous and implicit solution of both the stellar struc- 
ture equations and the chemical composition equations, etc. 
make the code very stable and easy to use. The current code 
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uses an equation of state that includes pressure ionization 
and Coulomb interaction, recent opacity tables derived from 
Rogers & Iglesias (1996) and Alexander & Ferguson (1994) 
(see Chen & Tout, 2007), nuclear reaction rates from Caugh- 
lan et al.((T985) and Caughlan & Fowler (1988), and neutrino 
loss rates from Itoh et al. (JTHHH HSHt • 

The ratio of mixing length to the local pressure scale- 
height a = l/Hp is set to 2, which fits to the Sun 
(|Pols et al. 1998). Convective overshooting is included by in- 
corporating a condition that mixing occurs in a region with 
V r > V a - S ov /(2.5 + 20£ + 16£ 2 ), where £ is the ratio of 
radiative pressure to gas pressure and <5 ov is a specified con- 
stant. 8 ov = 0.12 gives the best fit to the observed systems 
([Schroder ct al. 1997), which corresponds to an overshoot- 
ing of about 0.25-Hp. RLOF is included from the boundary 
condition 



CMaxfO, (^a 

r lobc 



-I) 3 



(2) 



dm _ 
dt 

where r sta r and ri b e are the radii of the mass donor and 
its Roche lobe, respectively, dm/dt gives the mass loss rate 
of the primary, and C is a constant. In our study, we set 
C = 5OOAf0yr _1 , with which RLOF proceeds steadily and 
the lobe-filling star overfills its Roche lobe as necessary but 
never overfills its lobe by much, typically (r s tar/n bc — 1) < 
0.001. 

If mass transfer is non-conservative, i.e. j3 is less than 
1.0, some matter will be lost from the system, taking away 
some angular momentum. In our study, the accretor is a 
main-sequence star, and it is much more compact than 
a giant. So we assume that the lost mass carries away the 
same specific angular momentum as pertains to that of the 
accretor, similar to the case of a compact component. The 
study of Beer et al. (|200Tp shows that this assumption is 
appropriate for main-sequence components. 

We only follow the evolution of the primary (initially 
more massive component) in a binary system when we 
study the critical mass ratio, since the structure of the pri- 
mary is the main factor for the formation of CE in this 
phase. We adopted the solar metallicity (Z = 0.02) 
in our calculations. As we are mainly concerned with 
BSs which originate from mass transfer, only low- 
and intermediate-mass binaries are studied here. The ini- 
tial mass of the primary increases from 1 to 8Mq by step 
of about AlogM/M = 0.1 (i.e. Mu=1.00, 1.26, 1,58, 2.00, 
2.51, 3.20, 3.98, 4.95, 6.31 and 7.94M Q , respectively). When 
the mass donor evolves to the giant branch, it dramatically 
expands, and then its radius R may well represent the evo- 
lutionary phases. So we systematically vary the radius of the 
mass donor at the onset of RLOE[j and the mass of the com- 
panion star for each primary mass. If Afu < 2.00Mq, the 
primary has a degenerate He core and undergoes a He flash 
as the central He is ignited. Before the He flash, the star 
can approach a very high luminosity as well as a very large 

3 The radius of the primary at the onset of RLOF is 
actually its Roche lobe R C ii, since the primary is just 
filling of its Roche lobe at that time. If the mass ratio 
q is given, the corresponding initial orbital separation A 
can be calculated from R crl /A = 0A9q 2 / 3 /[0.6q 2 / 3 + ln(l + 
q 1 ' 3 )] (Eggleton, 1983). For convenience, we only show 
the radius of the primary at the onset of RLOF in this 
paper. 




Figure 2. Evolutionary tracks for single stars from 1 to 8Mq. For 
primaries in binary systems, the dots (for FGB) and asterisks (for 
AGB) indicate their positions at the onset of RLOF we studied 
in the paper. Corresponding stellar radii are presented in Table 
[T] The AGB cases for Mi < 2Mq are not plotted to ensure the 
distinction of the figure. 



radius. So we choose an increase of the radius AlogR/R® by 
steps of 0.2 from the base of the first giant branch (FGB) 
to the tip of the FGB to study the critical mass ratio. When 
Mil > 2.OOM0, He burns quietly in the center and the ra- 
dius difference of the primary from the base to the tip of 
FGB is not very large. We therefore set Alogi?/ Rq = 0.1 
for this case. 

If a star has not filled its Roche lobe before central He 
burning, it will not fill its Roche lobe during He burning, 
because the star will contract after He ignition. We there- 
fore have not considered the He-burning phase in the stud- 
ies. After the exhaustion of its central He, the star expands 
again and its radius may become larger than that at the tip 
of the FGB. If Mu > 2M , we continue the study on the 
AGB when the stellar radius equals to that at the tip of the 
FGB, and increase the radius by steps of AlogR/R® = 0.1. 
Fig.[2]shows the positions of the primaries as RLOF begins, 
corresponding stellar radii are presented in Table 1. 

If Ma < 2Mq , the code breaks down when the He flash 
occurs. To investigate the case on the AGB for a mass donor 
with Mi < 2.OM0, we constructed some stellar models after 
the He flash. Considering that some of the mass donors 
on the AGB may have evolved from some initially more 
massive stars Q we choose the point at the minimum stellar 
radius at the tip of the FGB for stars with Mu > 2Mq 
(i.e. logR/RQ = 1.41) to start our studies and increase the 
radius by steps of A\ogR/RQ — 0.1. The stellar radii and 
corresponding q c are listed in Table 2. 
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Table 1. Critical mass ratio for stable RLOF. Afio is ZAMS mass for the primary; logfj is the radius of the primary at the onset 
of RLOF, q c is the critical mass ratio and /3 is the fraction of the mass lost from the primary accreted by the secondary. 



Afio 



logf? 



9c 
,8=0.0 



9c 

/3=0.5 



9c 

13=1.0 



Mio logR 



Qc 
3=0.0 



9c 
3=0.5 



9c 
3=1.0 



Mm 



log/? 



9c 
,9=0.0 



9c 
/3=0.5 



9c 

,3 = 1.0 



1.00 0.3048 


1.3943 


1.1927 


1.0607 


2.51 


0.9785 


1.4614 


1.2539 


1.0886 2.51 


1.6653 


1.2490 


1.0574 


0.8920 


0.5133 


1.2477 


1.0546 


0.8933 




1.0867 


1.2757 


1.0842 


0.9279 


1.7653 


1.2183 


1.0303 


0.8670 


0.7133 


1.2416 


1.0530 


0.8918 




1.1867 


1.2366 


1.0473 


0.8894 


1.8653 


1.1900 


1.0032 


0.8415 


0.9133 


1.2615 


1.0709 


0.9066 




1.2867 


1.2138 


1.0263 


0.8655 


1.9653 


1.1535 


0.9674 


0.8079 


1.1133 


1.2642 


1.0719 


0.9040 




1.3867 


1.1945 


1.0087 


0.8492 


2.0653 


1.1176 


0.9325 


0.7774 


1.3133 


1.2545 


1.0595 


0.8902 




1.4867 


1.1759 


0.9999 


0.8341 


2.1653 


1.0847 


0.9002 


0.7471 


1.5133 


1.2228 


1.0312 


0.8660 












2.2653 


1.0639 


0.8776 


0.7259 


1.7133 


1.1901 


0.9968 


0.8324 


3.20 


1.2151 


1.5982 


1.3970 


1.2560 


2.3653 


1.0512 


0.8652 


0.7168 


1.9133 


1.1403 


0.9561 


0.7958 




1.3236 


1.3360 


1.1371 


0.9678 


2.4653 


1.0340 


0.8505 


0.7053 


2.1133 


1.1131 


0.9257 


0.7696 




1.4236 


1.2672 


1.0765 


0.9128 




















1.5236 


1.2200 


1.0338 


0.8738 3.20 


1.7076 


1.2668 


1.1075 


0.9235 


1.26 0.4110 


1.3598 


1.1588 


1.0411 




1.6236 


1.1852 


1.0029 


0.8449 


1.8076 


1.1989 


1.0138 


0.8534 


0.6194 


1.2172 


1.0215 


0.8658 












1.9076 


1.1650 


0.9814 


0.8237 


0.8194 


1.2226 


1.0272 


0.8685 


3.98 


1.4407 


1.7427 


1.5413 


1.3747 


2.0076 


1.1236 


0.9417 


0.7880 


1.0194 


1.2341 


1.0433 


0.8783 




1.5492 


1.3730 


1.1754 


1.0029 


2.1076 


1.0871 


0.9059 


0.7549 


1.2194 


1.2266 


1.0360 


0.8704 




1.6492 


1.2796 


1.0892 


0.9238 










1.4194 


1.2089 


1.0174 


0.8532 




1.7492 


1.2172 


1.0308 


0.8713 3.98 


1.8793 


1.2251 


1.0347 


0.8737 


1.6194 


1.1704 


0.9814 


0.8208 




1.8492 


1.1704 


0.9880 


0.8305 


1.9793 


1.1827 


0.9992 


0.8413 


1.8194 


1.1289 


0.9399 


0.7835 












2.0793 


1.1399 


0.9581 


0.8043 


2.0194 


1.0871 


0.8983 


0.7471 


4.94 


1.6671 


1.7617 


1.5597 


1.3637 


2.1793 


1.1023 


0.9205 


0.7683 












1.7756 


1.3747 


1.1770 


1.0033 


2.2793 


1.0647 


0.8804 


0.7356 


1.58 0.5715 


1.3634 


1.1602 


1.0055 




1.8756 


1.2719 


1.0823 


0.9161 


2.3793 


1.0394 


0.8570 


0.7116 


0.7800 


1.2090 


1.0225 


0.8649 




1.9756 


1.1952 


1.0116 


0.8512 










0.9800 


1.2047 


1.0098 


0.8535 










4.94 


2.0583 


1.2105 


1.0198 


0.8628 


1.1800 


1.2094 


1.0202 


0.8583 


6.31 


1.9044 


1.7582 


1.5419 


1.3578 


2.1583 


1.1619 


0.9792 


0.8249 


1.3800 


1.1908 


1.0053 


0.8441 




2.0129 


1.3547 


1.1576 


0.9844 


2.2583 


1.1201 


0.9392 


0.7870 


1.5800 


1.1629 


0.9760 


0.8162 




2.1129 


1.2563 


1.0689 


0.9034 


2.3583 


1.0823 


0.9013 


0.7524 


1.7800 


1.1174 


0.9310 


0.7764 




2.2129 


1.1548 


0.9740 


0.8178 










1.9800 


1.0771 


0.8886 


0.7393 










6.31 


2.2613 


1.1864 


0.9971 


0.8392 










7.94 


2.1091 


1.6967 


1.4669 


1.2748 


2.3613 


1.1438 


0.9622 


0.8090 


1.91 0.7239 


1.3905 


1.1847 


1.0165 




2.2176 


1.2971 


1.1004 


0.9288 


2.4613 


1.0985 


0.9177 


0.7686 


0.9324 


1.2140 


1.0269 


0.8695 




2.3176 


1.2306 


1.0457 


0.8804 


2.5613 


1.0651 


0.8841 


0.7357 


1.1324 


1.1971 


1.0078 


0.8504 




2.4176 


1.1247 


0.9476 


0.7924 










1.3324 


1.1816 


0.9977 


0.8369 










7.94 


2.4438 


1.1624 


0.9764 


0.8201 


1.5324 


1.1652 


0.9803 


0.8217 




For 


AGB 


case 




2.5438 


1.1235 


0.9418 


0.7913 


1.7324 


1.1223 


0.9387 


0.7838 


2.00 


1.4175 
1.5175 


1.2764 
1.2661 


1.0819 
1.0730 


0.9130 
0.9063 


2.6438 


1.1019 


0.9203 


0.7700 


2.00 0.7672 


1.3710 


1.1779 


1.0065 




1.6175 


1.2388 


1.0474 


0.8805 










0.8759 


1.2447 


1.0530 


0.9035 




1.7175 


1.2068 


1.0171 


0.8526 










0.9757 


1.2261 


1.0298 


0.8745 




1.8175 


1.1698 


0.9816 


0.8210 










1.0757 


1.2039 


1.0157 


0.8614 




1.9175 


1.1329 


0.9470 


0.7890 










1.1757 


1.1980 


1.0093 


0.8521 




2.0175 


1.0971 


0.9090 


0.7546 










1.2757 


1.1890 


1.0030 


0.8446 




2.1175 


1.0901 


0.9045 


0.7462 










1.3757 


1.1786 


0.9936 


0.8355 




2.2175 


1.0734 


0.8844 


0.7313 











3 THE CRITICAL MASS RATIO FOR A 
GIANT MASS DONOR 

As mentioned in section 1, dynamically unstable RLOF of- 
ten occurs when the mass donor is a giant at the onset of 
mass transfer. Fig. [3] shows the mass transfer rate vs the 
mass of the primary for a typical binary undergoing dynam- 
ical unstable mass transfer. The primary is l.OOMg with a 
core mass M c = O.344M at the onset of RLOF, and /3 — 



4 For example, a star with Mij > 2Mq becomes less than 2Mq 
because of stellar wind, then the star has not passed through a 
phase with a high radius before He burning. 



during the RLOF. From the figure we see that mass transfer 
initially occurs on a thermal time-scale. The mass transfer 
rate, M, rises quickly to a level of 10 _5 AfQyr _1 at first and 
continues to grow more slowly later. After the primary 
has lost about O.CMM© (point A), the system encounters 
dynamical instability. The behaviour here is similar to that 
of radiative mass donors with a mass ratio greater than 3 
(Han & Podsiadlowski 2006), but it cannot be explained by 
the evolution of the entropy profile as in the case of a mass 
donor with a radiative envelope. It is related to the evolution 
of binary parameters, which strongly depend on the angu- 
lar momentum loss associated with the mass loss from the 
system. Furthermore, in comparison to the delayed dynami- 
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Table 2. Critical mass ratio for stable RLOF for AGB mass donors with mass less than 2Mq. M\ is the mass of the primary, 
logR is the radius of the primary at the onset of RLOF, q c is the critical mass ratio and (3 is the fraction of the mass lost from the 
primary accreted by the secondary. 



logR Mi = 1.00M Q 



Mi = 1.26M, 







Mi = 1.6OM 



B., 



(3 = 0. 



13 = 0.5 (3 = 1. 



(3 = 0. 



13 = 0.5 (3=1. 



(3 = 0. 



(3 = 0.5 (3 = 1. 



1.4096 


2.0406 


1.7919 


1.5816 


1.7504 


1.4589 


1.2640 


1.5018 


1.3020 


1.1192 


1.5096 


1.9040 


1.6705 


1.4550 


1.5912 


1.3709 


1.2294 


1.4360 


1.2290 


1.0533 


1.6096 


1.8061 


1.5761 


1.3728 


1.5302 


1.3154 


1.1481 


1.3954 


1.1904 


1.0157 


1.7096 


1.7154 


1.4891 


1.3096 


1.4791 


1.2676 


1.0833 


1.3617 


1.1608 


0.9834 


1.8096 


1.6223 


1.4045 


1.2144 


1.4142 


1.2078 


1.0256 


1.3059 


1.1100 


0.9362 


1.9096 
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Figure 3. Mass transfer rate vs the primary mass for a binary 
of Mn = I.OOMq. Mass transfer efficiency (3 = 0.0 and the core 
mass is 0.344Mq. Point A gives the place at which mass transfer 
dramatically increases because of dynamical instability. 



cal instability in the radiative mass donor, the primary with 
a convective envelope here loses much less mass during the 
delayed time and probably has little influence on the final 
outcome. 

The results of our calculations are summarized in Tables 
[T] and [2] For each ZAMS, we present the radius of the pri- 
mary logR/R® at the onset of RLOF and the critical mass 
ratio q c for dynamically stable RLOF from different mass 
transfer efficiency (3. The results demonstrate that RLOF is 
probably dynamically stable even if the mass donor is more 
massive than the secondary for (3 < 1. In particular, when 
the mass donor is less than 2Mq and fills its Roche lobe on 
the AGB, RLOF may be dynamically stable even for (3 = 1. 
Since the core is not very degenerate and the envelope is not 
fully convective when the star is at or near the base of the 
FGB, q c is much larger at the base of the FGB than in the 
following evolutionary phases for each /3. 

Figs.|4]and[5]present q c vs logR/RQ— Ao for low-mass 
and intermediate-mass binaries, respectively, when the 
mass donors are on the FGB. Here Ao is the log of the 
radius of the primary at the base of the FGB. We 



include the result of Mu = 1.9M© in Fig. 0] (see also Table 
[l]). The results for the mass donor on the AGB are showed in 
Figs.[6]and[7] In all of the figures, we see that a more evolved 
star (with a larger stellar radius at the onset of RLOF) is 
less stable for RLOF, especially when the mass donor is 
on the AGB at the onset of mass transfer. The reasons are 
as follows. For more evolved binaries, the evolutionary time- 
scale is shorter and hence the mass transfer rate higher than 
those for the less-evolved ones. However q c is non-monotonic 
initially in Figf3]because the core is not very degenerate and 
the envelope is not yet full convective in this phase. 

The cases at the last point on the FGB and at the 
first point on the AGB give a different tendency as de- 
scribed above, since the two points are on different giant 
branches, and the stars have different cores, i.e. a non- 
degenerate He core on the FGB while a possibly de- 
generate CO core on the AGB. Meanwhile, the degree 
of convection is also different at the two points. 

As shown in Fig0 there is a threshold for q c when the 
initial primary mass Mu - 3.98M Q . When Mr, < 3.98M , 
g c systematically increases with the initial primary's mass 
while it is opposite when Mu > 3.98A/q. This phenomena 
is likely revelant to the degenerate degree of the core. We 
examined the degeneracy parameter ip in the calculations 
and found that tp becomes negative just at Mu = 3.98Mq. 

We know that, when central He is exhausted in a star, 
there exists a He-burning shell and a H-burning shell. The 
H-burning shell extinguishes first when it burns outwards, 
and then the He-burning shell also extinguishes. The star 
begins to contract and the temperature increases until H is 
ignited. With the increase of temperature, He is also ignited 
but in a very thin shell. The He ignition in a very thin shell 
is unstable, and leads to a dramatic expanding of the star. 
After the expansion, He-burning becomes stable and the star 
then has two stable burning shells again, and so on. This 
behaviour is known as the thermal pulse. The thermal pulse 
may result in a sharp increase of mass transfer rate and leads 
dynamical instability to be more likely to occur such as 
in Fig. [8] In principle, if the mass donor has two burning 
shells, the thermal pulse will likely affect the final results. 
For the cases we studied, however, the influences seem small 
and can be neglected. 
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Figure 4. Critical mass ratio for low-mass binaries of FGB+MS, 
i.e. Mu = 1.00, 1.26, 1.60 and 1.9OM . The lines are from the 
fitting formulae of equations (4) to (6). 

To simply use this results, we fitted q c from our calcula- 
tions. The influence from (3 is nearly linear and can be fitted 

as 

<7c = (fco — co/3, (3) 

For the FGB, 

q c0 = ci + c 2 QogR/Ro - AG) + c 3 (logR/R Q -A0) 2 , (4) 

when Mu < 2.00M© and logR/R e - A0 > 0.2. In other 
cases, 



(•») 



g c o = c 4 + c 5 {lo S R/R Q - A0) 1/3 . 

Here 

ci = ci,i + Ci, 2 Mi + c i>3 Mi . (6) 

The values of Cij and Co are listed in Table [3] for different 
cases. 

For the case of the AGB, 

q c0 = ci + c 2 (logR/Rs - A0) + c 3 (logR/R @ - AO) 2 . (7) 

The coefficients Ci and Co are listed in Table [4] One 
may send a request to xuefeichen717@hotmail.com for the 




1.5 
logR/R 3 -A 

Figure 6. Critical mass ratio for low-mass binaries when the 
primary is on the AGB. The lines are from the fitting formulae 
of equations (6) and (7). 



FORTRAN code of all the formulae used in this pa- 
per. 

To compare with the results of the polytropic model, 
we present the dependence of q c on the core mass M c as 
well as on the core mass fraction f c for low-mass binaries in 
Fig. [9] ignoring the points at the base of the FGB for each 
mass. Here the core mass is defined as the mass within the 
hydrogen mass fraction less than 0.1 q From the figure, we 



5 It is a little bit difficult to determine the core in a real star 
as described in the polytropic model. However, the core mass 
defined here is a close approximation to that denned in 
the polytropic model. 
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Figure 7. Critical mass ratio for intermediate- mass binaries when 
the primary is on the AGB. The lines are from the fitting formulae 
of equations (6) and (7). 
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Figure 9. The dependencies of critical mass ratio q c on the core 
mass fraction f c and on the core mass M c for low-mass binaries, 
i.e the initial primary mass Mu = 1.00, 1.26, 1.6 and 1.9Mq, 
respectively. The triangles, dots and asterisks are for mass transfer 
efficiency f3 = 0., 0.5 and 1.0, respectively. The squares and pluses 
are from Han et al (2002), where stellar wind is included before 
RLOF and /3 = 0. The pluses are for binaries with a primary mass 
equal to O.8M , and the squares for M i; > l.OOM . The solid 
lines in the bottom panel are from equations (8) and (9). 



Table 3. The coefficients c; j and cq for various cases when 
the mass donor is on the FGB (from equations(7) to (10)). 
Here r = logR/R e - A0. 



cases 


i 


i = i 


J = 2 


i = 3 


CO 


Mi < 2.00M Q 

& 

r > 0.2 


1 
2 
3 


1.4418 
-0.08252 
0.08967 


-0.3477 
0.4180 
-0.3129 


0.1180 
-0.1873 
0.1079 


0.37 


Mi < 2.OOM 
& r < 0.2 


4 
5 


1.7141 
-0.3978 


-0.4919 
0.2620 


0.1695 
-0.1101 


0.37 


Mi > 2.OOM 
& Mi < 3.98M 


4 
5 


0.9989 
0.2740 


0.1755 
-0.2353 


0.003139 
-0.007348 


0.35 


Mi > 3.98M© 


4 
5 


1.3840 
0.05596 


0.1389 
-0.2602 


-0.0124 
0.02045 


0.37 



Figure 8. Mass transfer rate as well as H-burning luminosity vs 
the primary mass for a binary when the mass donor is on the 
AGB. 



see a clear relation between q c and M c , but not between q c 
and f c . q c may be fitted as follows: 



q c — q c o - 0.35/3 

where 

q c0 = 1.142 + 1.081M C - 2.852M C 2 



(8) 
(9) 



The studies above have not included the influences from 
the stellar wind. Han et al. l(2002i included the stellar wind 
prior to RLOF when they studied the q c for low-mass bina- 
ries on the FGB and j3 — 0. We have included their results 
in Fig. |9] The pluses are for binaries with Mn = 0.8M(^ 



° We have not studied the case for Mn = 0.8Mq in this pa- 
per, because the time-scale is too long for a star with this mass 
to evolve from zero-age main sequence to giant branch. Mean- 
while, it might be more likely that low-mass binaries (i.e 
the initial primary less than IMq) contribute to blue 
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Table 4. The coefficients Cij and eo for various cases when 
the mass donor is on the AGB. 



cases 


i 


J' = l 


i = 2 


J = 3 


CO 


Mi < 2.51M 


1 
2 
3 


1.0462 
-0.6179 
0.7346 


0.3115 
0.5193 
-0.7568 


-0.03702 
-0.2027 
0.2074 


0.34 


Mi > 2.51M Q 

& 
Mi < 3.98M 


1 
2 
3 


-0.5270 
7.1432 
-4.9258 


1.4289 
-4.9954 
3.2653 


-0.2325 
0.7625 
-0.4966 


0.36 


Mi > 3.98M Q 


1 
2 
3 


1.4681 
0.1092 
-0.1971 


0.01922 
-0.2844 
0.1226 


-0.004501 

0.02368 

-0.006051 


0.36 


Mi < 2.00M Q 


1 
2 
3 


0.02442 
9.6761 
-4.5199 


7.0840 

-19.0892 

7.8189 


-3.6961 
7.9462 
-3.1087 


0.39 



and the squares for Mu > LOOM© (see Table 3 in their pa- 
per). We see that all the squares are well along the fitting 
line of f3 = 0, indicating that stellar wind has little influence 
on q c . However, the mass ratio will decrease due to stellar 
wind, and RLOF will be more stable, e.g. the mass ratio 
possibly becomes less than q c at the onset of RLOF and the 
binary will undergo a stable mass transfer. For the case of 
Mu - O.8M , q c from Han et al. (20021) is obviously larger 
than the fitting value, which means that RLOF is more sta- 
ble in binaries with this primary mass. Since the binaries 
with Mu = 0.8M@ have a larger f c in comparison to those 
with Mu > 1.0M© Q the high q c for M u = 0.8M Q here 
likely indicates some contributions from / c , as is possibly 
concealed by that from M c for binaries with Mu > LOOM© 
(see the discussion in section 6). 



4 EXAMPLES OF THREE BINARIES 
RESULTING IN BSS 

In this section, we present the detailed evolutions for three 
binaries, which are selected mainly for the reason that 
we are interested in BSs in the old open cluster M67, where 
many BSs, including several long-period BSs, have been ob- 
served. M67 has a metallicity similar to that of the 
Sun (|Hobbs fc Thorburn 199UIFriel fc Janes 1993|) . 
There are some researches on the age of M67. It 
may range from 3.2+0.4 Gyr (jBonatto fc Bica 20031) 
to 6.0 Gyr Panes fc Phelps 199~4l ). The study of 
VandenBerg & Stetson (2004) derived an age of 
4.0 Gyr. In the N-body model of this cluster 
(Hurley et al. 2005), the authors investigated the 
behaviour around 4 Gyr. So we also consider that 
its age is about 4 Gyr, indicating that the mass of 



stragglers via coalescence induced by angular momen- 
tum loss (Chen &: Han 2008), while not from the mass 
transfer between a giant and a MS. 

7 In the study of Han et al. (200.2), the mass donor has a similar 
core mass but a different stellar mass at the onset of RLOF. 
Therefore, a high stellar mass means a smaller core mass fraction 
fc- See Table 3 in their paper for details. 



the turnoff M to is about 1.2 to 1.3Mq, which is the 
main factor of the cluster we will consider as we 
choose the binary samples. Basic requirements for the 
binaries which may contribute to BSs in a cluster via a giant 
transferring mass to a MS companion are for the masses of 
both components, i.e. the primary should be more massive 
than the turnoff to ensure that it has left the main sequence 
and the secondary should be less than the turnoff if it is still 
on the main sequence at the cluster age. The fact that the 
secondary should accrete certain material before it becomes 
a BS will give a further constraint on the binary parameters. 
Obviously, j3 = is ruled out since no material is accreted by 
the secondary in this case. In the following of this section, 
we will give three examples and present the evolutionary 
results in details. 

Example 1: RLOF is dynamically stable for initial pa- 
rameters. As shown in section 3, q c decreases with the mass 
transfer efficiency j3. This means that, in order to ensure 
that /3 is large enough to increase the secondary's mass to 
be larger than the turn-off, the mass ratio q should be as 
small as possible. However, q should be larger than unity to 
confirm that mass transfer occurs between a giant star and a 
MS companion. We therefore choose a binary of 1.3 + 1.2Mq 
(q = 1.08). Meanwhile, since the orbital period increases 
with the core mass (or the radius) of the mass donor at 
the onset of RLOF, a large core mass (or stellar radius) 
is necessary for long orbital period BSs. However, when 
M c > 0.4M Q , q c is less than 1.0 even for (3 — 0. We there- 
fore set M c = 0.35M© (logR/R Q = 1.7503) at the onset of 
RLOF. From equations (8) and (9), (3 = 0.16 for q c = 1.08. 
We then set j3 = 0.1 in the calculation. 

Example 2: RLOF is dynamically unstable for initial 
parameters, but it will be stable after some mass loss of 
the primary by stellar wind prior to RLOF. Stellar wind is 
included only on the giant branch by the mode of Reimers' 
lfT975l : 



AU 



4 x 10~ 13 rjRL/M, 



(10) 



where r\ is a dimensionless factor. The same binary as 
example 1 was examined here but with j3 — 0.5. In gen- 
eral, rj is 1/4 (|Renzini 198H llben fc Renzini 19831 
ICarraro et al. 1996[) . However it is expected by 
many authors that giant stars have much higher 
mass loss rates than that of Reimers, especially 
when the stars are far away from the base of giant 
branch (e.g. Bloecker 1995). Furthermore, tidal in- 
teraction between the two components of the binar 
also increases the mass loss rate. So we set 77 = 2.5 
in the calculation. When RLOF begins (M c = 0.355M Q ), 
Mi = l.llM e (q = Mi /M a = 0.93). We switch the stellar 
wind off once the mass transfer rate exceeds the value given 
by eq(10) and have not included it after RLOF. From eqs(8) 
and (9), we get q c = 0.99 for (3 = 0.5 and M c = 0.35M Q . So 
RLOF is dynamically stable as it begins. 

Example 3: RLOF is dynamically unstable for initial 
parameters, but it becomes stable after the initially dy- 
namically unstable mass transfer, if we assume that CE 



8 We also examined the cases for r\ = 1/4 and 1, and found that 
both of them are too small to strip enough mass away to lead to 
q < q c as RLOF begins. 
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Figure 10. Mass loss rate vs the mass of the primary for the three 
binaries in the text. The system parameters are presented in each 
panel. Stellar wind is not included in panels (a) and (c), while 
it is included at the giant branch in panel (b) with the mode of 
Reimer (r] = 2.5). During the initially dynamically unstable phase 
in panel (c), we set the highest mass loss rate of the primary to 
be 1 X lO _4 M0yr _1 (see the text for details). 



has not developed during this phase according to the study 
of Beer et al. (12007)) . Beer et al. I|2U07|) argued that a 
wide range of systems avoid CE phase because mass trans- 
fer is super-Eddington even for non-compact companions. 
The accretion energy released in the rapid mass-transfer 
phase strips away a large fraction of the giant's envelope, 
reducing the tendency to dynamical instability and merg- 
ing. More constraints are necessary to determine whether 
the CE is formed or not in this case, but we have little 
knowledge about it. For simplicity, we assume that CE has 
not developed in the binary we study here, i.e. in a binary 
with 1.6 + I.IMq. The code cannot work for the initial dy- 
namically unstable mass transfer because of the high mass 
transfer rate M. We therefore artificially limit the highest 
M to be 1 x 10~ 4 Moyr~ 1 , and the accretion rate of the 
secondary M a is equal to M as M < 1 X lO~ 5 M yr _1 , 
and equal to 1 x KT 5 Af yr" 1 as M > 1 x lO~ 5 M yr _1 . 
From this assumption, /3 ranges from 1.0 to 0.1, then to 1.0 
again during the whole RLOF period. The maximum ac- 
cretion rate of the MS companion here, i.e. lx 10 M yr _1 , 
is comparable with the typical value of symbiotic stars 
HScott fc Webbink 1984|) . 

Figure [10] presents the mass loss rate to the total mass 
of the primary for the three binaries. In these figures we see 
that, after the initial rapid mass loss, the primary loses its 
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Figure 11. The evolutionary tracks for the secondaries of the 
three binaries we examined in the text. 



material at a rate of about 10~ to 10 _ M@yr~ . Generally 
this slow mass transfer occurs after the inversion of the 
mass ratio, i.e. the primary is less massive than the sec- 
ondary. During the slow mass transfer phase, the radius of 
the primary again increases but with a much lower level 
in comparison to that in the rapid mass loss phase. At the 
same time, the Roche lobe radius of the primary increases 
with mass loss, since the orbital period becomes longer and 
longer after the inversion of q. So if CE has not developed 
in the dynamically unstable phase in a binary such as ex- 
ample 3, it is very likely that CE will never develop in the 
binary. The material lost from the primary in this phase will 
be lost from the system in a way similar to stellar wind, 
and RLOF eventually terminates after most of the primary's 
envelope is lost. Note that the mass loss of the system here 
is different from the ejection of CE. The former needs no or- 
bital energy but the latter needs some. So binaries avoiding 
CE have long orbital periods. 

The behaviours of the secondaries of the three binaries 
are presented in Fig. 1111 Since the secondaries cannot ac- 
crete the mass from the primary in a very short time, they 
have left thermal equilibrium and evolve in a way similar to 
pre-MS stars during the initial phase instead of going up- 
ward along the main sequence. Due to the long orbital peri- 
ods, the secondaries do not overfill their Roche lobes. They 
finally become normal main-sequence stars with a higher He 
fraction in the envelope when new thermal equilibrium is es- 
tablished, and evolve similarly to normal single stars with 
those masses. These rejuvenated stars have much longer 
timescales (including the phase of the secondaries with lower 
masses before accretion) on the main sequence and have the 
possibility to be BSs. 

The main characteristics of the three binaries are listed 
in Table from which we may obtain some clues on their 
parent binaries and evolutionary histories. For example, al- 
though the secondary in example 2 has a mass similar to 
that in example 3 after accretion, the latter is much bluer 
(Fig. Ill[) when both of them return to thermal equilibrium 
since the latter is much less evolved at the onset of RLOF. 
The two secondaries have different orbital periods, different 
compact components and different lifetimes etc., and appear 
in different age of the cluster, as shown in Table [5] 
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for simplicity, we only consider a constant mass ratio distri- 
bution (jMazeh et al,1992[) . 




-0.5 



(B-V) 

Figure 12. The colour and magnitude evolutions for the indi- 
vidual components of the binary 1.6 + 1.1-Mq as well as for the 
whole system. The early evolutions of both components prior to 
RLOF are not plotted in the figure. The circle and dotted lines 
are for the primary and the secondary, respectively, and the solid 
line is the combination of the two components. The numbers 
are the time sequences of the evolutions. 



As an example, we give a colour-magnitude diagram of 
both of the components as well as the binary system for 
example 3 in Fig |12l The combined colour and magnitude 
of the binary result from add-and-subtract calculations of 
the luminosity of both components (see Appendix A). The 
evolutions of the components prior to RLOF have no differ- 
ences from a standard stellar evolution and are not plotted 
in the figure. During the whole mass transfer phase the sys- 
tem resembles the characteristics of the primary because the 
primary is much more luminous than the secondary, but the 
timescale is very short. Once the primary becomes cooler, 
i. e. less luminous than its component, the binary is hardly 
distinguishable from the main-sequence component. It is 
very difficult to find the clues from the compact component, 
even from spectral observations because of the long orbital 
period. 



5 MONTE CARLO SIMULATION 

To investigate BSs from mass transfer between giants and 
MS companions, we performed a Monte Carlo simulation 
for a sample of 10 6 binaries (very wide binaries are actually 
single stars). A single starburst is assumed in the simulation, 
i.e. all the stars have the same age and metallicity (Z — 
0.02). The initial mass function (IMF) of the primary, the 
initial mass ratio distribution and the distribution of initial 
orbital separation are as follows: 

i) the IMF of Miller & Scalo (fT979J) is used and the 
primary mass is generated from the formula of Eggleton, 
Fitchett & Tout JTSB5JI: 



Mi = 



(l-X)°- 75 +0.032(l-X) 1 /4 



(11) 



where A is a random number uniformly distributed between 
and 1. The mass ranges from 0.8 to WOMq. 

ii) the mass ratio distribution is quite controversial and, 



n{q) = l,0<g' <1 



(12) 



where q' = l/q - M 2 /Mx. 

iii) We assume that all stars are members of binary 
systems and the distribution of separations is constant in 
logo (where a is separation). 



an(a) = 



,(a/ao) r 



a < ao 

ao < a < a\ 



(13) 



where a = 0.070, a = 10.R©,ai = 5.75 x 10 6 i? Q = 0.13pc 
and m — 1.2. This distribution gives an equal number 
of wide binary systems per logarithmic interval and 50 
per cent of systems with orbital periods less than 100 yr 
Pl (|Han et al. 1995[> . 

The rapid binary evolution code developed by Hurley 
et al. (2000; 2002) is employed here. In addition to all as- 
pects of single-star evolution, this code includes many fea- 
tures of binaries, i.e. mass transfer, mass accretion, common- 
envelope evolution, collisions, supernova kicks and angular 
momentum loss mechanisms etc.. In particular, circulariza- 
tion and synchronization of the orbit by tidal interaction are 
calculated for convective, radiative and degenerate damping 
mechanisms. As a comparison, we adopt three criteria ob- 
tained in different ways, i.e. from a polytropic model (equa- 
tion 1), from the paper of Hurley et al. (2002) and from this 
paper. 

In the paper of Hurley et al. |2002 ) , 
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[1.67 - x + 2(M C /M) 5 ]/2.13, 



(14) 



where M c and M are the core mass and the total mass of 
the donor, respectively, x is the exponent of the mass-radius 
relation at constant luminosity for giant stars and equals to 
0.3. This criterion comes from the assumption that the adia- 
batic mass-radius exponent of a giant star £ a( j = d\nR/d\nM 
equals to the Roche lobe mass-radius exponent £l, where 
£l ~ 2.13q — 1.67 for conservative mass transfer (Tout et 
al. 1997) and ( ad fa -x + 2(M C /M) 5 , which is fitted from 
detailed stellar models. 

From the sample, we obtain 303291 binaries which be- 
gin RLOF when the primary is a giant or more evolved star 
via Hurley's rapid binary evolution code. In the following 
section we will give the consequences from different cri- 
teria of the dynamical instability for mass transfer from a 
giant to a MS companion. 



5.1 The results from Hurley's criterion 

Among the 303291 binaries undergoing mass transfer from 
a giant or more evolved star to its companion, most of them 



9 It is an assumed distribution inferred from observations of spec- 
stropic binaries (see a series of papers of Griffin, R. F in The 
Observatory). Indeed, there are many distributions used in the 
literature, but they are all flat for wide binaries, leading to sim- 
ilar binary population results. The important thing for the cur- 
rently adopted distribution is that it implies 50 per cent of stellar 
systems with orbital periods less than 100 yr. One can simply 
multiply the results with a coefficient if another percentage is 
assumed. 
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Table 5. The characteristics for three different evolutionary examples in the text. 77 = the factor of Reimer's wind, 
M c = the core mass at the onset of RLOF, /3 = the mass fraction of the lost matter of the primary accreted by 
the secondary, i^LOF = the age at the onset of RLOF, £ms = the age of the secondary when it terminates its main 
sequence P e , M\ e , M2 e are the orbital period, the mass of the primary and the mass of the secondary at the end of 
RLOF. 
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example 1 1.3+1.2 0.0 0.356 0.1 803 4.71 5.45 0.46 1.28 

example 2 1.3+1.2 2.5 0.356 0.5 722 4.71 5.45 0.44 1.53 

example 3 1.6+1.1 0.0 0.340 - 439 2.47 4.57 0.41 1.50 



eventually experience CE evolution, leaving short-period 
binaries or mergers of the two components if the CE cannot 
be stripped awavl 10 ! According to Hurley's criterion, only 
3746 binaries may avoid CE formation and probably show 
the characteristics of some strange objects, e.g. symbiotic 
stars, BSs etc. during or after RLOF. Here we are concerned 
with the outcomes of BSs. 

Most of the 3746 binaries which avoid CE formation be- 
gin RLOF when the primaries are on the giant branches, i.e. 
on the first giant branch(FGB) for 17 binaries, on the early 
asymptotic giant branch (EAGB, after central He burning 
but before the first thermal pulse) for 415 binaries, and on 
the thermally pulsing AGB (TPAGB) for 3295 binaries. The 
simulation shows that, 3139 binaries will pass through BS 
phase during their lives and 2208 of them have undergone 
dynamically stable RLOF before they become BSa 11 !. This 
means that RLOF is an important process to increase the 
secondary's mass to be larger than the turnoff of a cluster. 
All the 2208 binaries have long orbital periods (greater than 
1000 d) and their mass ratios are less than unity at the on- 
set of RLOF, indicating that the primaries have lost most 
of their envelopes prior to RLOF (e.g 2164 binaries have a 
mass donor being on TPAGB). As a consequence, the BSs 
from this way also have long orbital periods (greater than 
1600 d). 

At the age of 4 Gyr, we obtain 96 BSs. Fig. 1 131 presents 
the distribution of the orbital period as well as the mass 
ratio versus the mass for the 96 BSs. In the figure, we see 
that, though the mass ratio at the onset of BS phase is a 
little different, all the BSs have similar mass ratios (around 
0.4) when they leave the main sequence. We explain this 
as follows. BS formation from mass transfer between 
giants and MS companion has some constraints. For 
example, RLOF should be dynamically stable. For 
q c adopted here, the primaries should be close to or 
on the TPAGB before 4 Gyr. Meanwhile, the sec- 
ondary should be stay on the MS at 4 Gyr (the MS 
timescale of the secondary after RLOF is only in or- 



10 Because the primaries have left the main sequence, these merg- 
ers are different from those from two main-sequence components. 
They are giant stars and then have no contribution to single BSs. 

11 We included wind accretion in the code for the reason pre- 
sented in section 6, so some BSs can be formed without RLOF. 
The simulation here indicates that wind accretion is important 
for BS formation, as the secondaries in some binaries cannot be 
larger than the turnoff mass without wind accretion even 
after RLOF . 
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Figure 13. Mass ratio as well as orbital period versus the mass 
of the secondary for 96 binaries which fulfill ti < AGyr < t%, 
where ti and t^ are the ages at the onset and at the end of BSs, 
respectively. The circles are for the case at t\ and the pluses are 
for the case at *2- 



der of 10 8 yr). So the binaries which may produce 
BSs in this cluster are very similar initially, and the 
products after RLOF are also similar. The differ- 
ence of q at the onset of BSs comes from the fact 
that RLOF has not terminated in these binaries at 
that time. Although RLOF still exists in some BSs, which 
results in the increase of their masses, the orbital period do 
not show obvious changes during the whole BS phase. All 
the BSs are near the turnoff of M67, indicating that this evo- 
lution channel, i.e mass transfer from a giant star to a main 
sequence companion(MS), can only account for the BSs in 
this region. 

5.2 The results from the polytropic model's 
criterion 

Using the same binary sample, we examined the evolution- 
ary consequence from the criterion of the polytropic model 
(see equation (1) in section 1). As discussed in the paper 
of Hurley et al. (2002), the value of q c from the polytropic 
model is obviously larger than that of Hurley et al. (2002) 
with increasing core mass. For example, it is a factor of 2 
larger at M c /Mi = 0.6. So more binaries may avoid dynami- 
cally unstable RLOF in the criterion of the polytropic model 
than those in the criterion of Hurley et al. (2002). 

From equation 1, 8522 binaries avoid CE formation and 
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5919 pass through the BS phase during their lives. Among 
of the 5919 binaries, 4287 have undergone or are just un- 
dergoing dynamically stable RLOF before the secondaries 
become BSs. At the age of 4 Gyr, we obtained 175 long- 
orbital period BSs (greater than 1600 d) . The distribution 
of the orbital period as well as the mass ratio versus the 
mass for the 175 BSs are similar to those from the criterion 
of Hurley et al. (2002). 



5.3 The results from this paper's criterion 

Now we study the outcomes from the criterion in this pa- 
per. Different from the two criteria above, we only 
evolve binaries to the onset of RLOF using Hurley's 
code and obtain the following evolutions from as- 
sumptions below. 

When we determine whether mass transfer leads to a 
BS formation or not in a certain cluster, two characteristic 
parameters of the secondary, i.e. the mass M a and the MS 
lifetimes £ms, are critical. The former is relevant to the ini- 
tial mass of the secondary and the mass transfer efficiency 
f3, while the latter is determined by the central H fraction 
and the total mass of the secondary after accretion. Mean- 
while, the ages at the onset and termination of RLOF 
are also important. For example, RLOF should start before 
or at the cluster age Cluster to ensure that the secondary may 
accrete some matter, but it cannot stop much earlier than 
tciuster, or the secondary has likely left the main sequence. 
Therefore, only the binaries, which complete mass transfer 
between £ c i us tcr — £ms and t c iustor as well as have secondary's 
mass M a > Mto during or after accretion, will contribute to 
BSs in a cluster. In fact, the mass transfer process from a gi- 
ant star to its companion is on a very short time-scale (about 
10 6 yr) in comparison to the MS time of the secondary after 
accretion, and then can be ignored. The amount of mat- 
ter lost from the mass donor during RLOF and the mass 
transfer efficiency f3 finally determine how much matter is 
accreted by the secondary. From section 4, we see that the 
final mass of the mass donor after RLOF is a little larger 
(about O.IM©) than the core mass at the onset of RLOF. 
Meanwhile, from the studies of the initial-final mass rela- 
tion ( |Meng et al. 2007| , the final mass of a star after RLOF 
or the super wind is around the core mass at the onset of 
RLOF or before super wind. The maximum difference in- 
duced by thermal pulses and mass loss is less than O.4M0, 
so we choose a range of the final mass of the mass donor 
from M c to M c + O.4M0 to examine the characteristics of 
the secondaries and of the binary systems. 

For M67 (truster = 4 x 10 9 yr), we obtained different 
possible candidates from various f3 for 4ms = 1 X 10 9 yr 
and £ms =5x 10 8 yr, respectivehl 12 !. The final orbital pe- 
riod of the candidates is estimated based on the assumption 



12 The MS time-scales here are referred to Table 5— it is about 
8 X 10 s yr from the onset of RLOF to the termination of 
the secondary on the main sequence. From Fig. 17, we see 
that the maximum mass ratio at the onset of RLOF is less than 
1.2 for binaries resulting in BSs, which means that the secondary 
is only slightly less than the turn-off, and then the MS time-scale 
of the secondary after accretion will not be very long. Meanwhile, 
it is likely longer than 5 X 10 s yr from Table 5. 



Table 6. The BS numbers in M67 obtained from FGB stars trans- 
ferring matter to MS companions for various 0. The MS timescale 
of the secondary after RLOF, <ms i i s se t to 1 X 10 9 yr (from the 2nd 
to 5th colomuns) and 5 X 10 s yr (from the 6th to 9th colomuns), 
respectively. The final masses of the primary are M^ = M c , 
M 2 f = M c + 0.1Af Q ,Mf t = M c + 0.2M Q and Mf t = M c + O.4M , 
where M c is the core mass of the primary at the onset of RLOF. 



ii 



M\ f M\ f 



Mf { Mf { AfL 



Ml M* M* 



0.1 414 412 406 371 319 317 311 276 

0.2 169 169 169 169 150 150 150 150 

0.3 31 31 31 31 30 30 30 30 

0.4 22222222 



Table 7. The BS numbers in M67 obtained from EAGB stars 
transferring matter to MS companions for various 0. 
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that the lost matter carries off the specific angular momen- 
tum as the accretor. Tables [6] to [8] present the BS num- 
bers for different mass donors, i.e. the mass donors are FGB 
stars (FGB+MS), EAGB stars (EAGB+MS) and TPAGB 
stars (TPAGB+MS), respectively. The corresponding distri- 
butions of the orbital period versus the final mass for various 
/3 are shown in Figs. H4ltoll6l 

We have not obtained BSs from an FGB star transfer- 
ring matter to a MS companion in M67 from the criteria 
of Hurley et al. (2002) and of the polytropic model, since 
the mass donor has not lost much mass by wind prior to 
RLOF and the core mass is not very large. Both of these 
facts make the condition q < q c difficult to fulfill (see 
equations (1) and (14)). However, from the criterion of this 
paper, several BSs will be formed in M67 if /3 < 0.4. From 
Table [6] we see that j3 strongly affects the contribution to 
BSs from this evolutionary channel. A small /3 contributes 
more BSs because RLOF is stabilized more easily. The mass 
of the secondary in this case is generally large enough to be 
larger than the turnoff after accretion even as /3 < 0.1. 

The contribution to BSs from EAGB+MS, in compar- 
ison to that of TPAGB-f-MS, is very small, although the 
value of q c for EAGB+MS is likely larger than that of 
TPAGB+MS because of the small radii at the onset of 
RLOF. This is relevant to the mass loss by wind prior to 
RLOF during the two phases. The stellar wind is generally 
much less in EAGB than that in TPAGB. Our study indi- 
cates that the mass loss prior to RLOF during EAGB phase 
is usually not enough to make q < q c , while a lot of mat- 
ter in the envelope of the primary has left the system by 
wind prior to RLOF for binaries of TPAGB+MS. So for 
systems of TPAGB+MS, the mass transfer is eas- 
ily stabilized and leads to the formation of long-orbital 
period binaries after RLOF. 
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Table 8. The BS numbers in M67 obtained from TPAGB stars 
transferring matter to MS companions for various /3. 
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Figure 14. Orbital period versus mass for the BSs in M67 re- 
sulting from FGB stars transferring material to main-sequence 
companions. The final mass of the primary is assumed to be 
M c + 0.2AfQ, where M c is the core mass of the mass donor at 
the onset of RLOF. See the text for the details. 
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Figure 15. Similar to Fig. 1141 but the mass donors are EAGB 
stars, i.e. after central He burning but before the first thermal 
pulse. 
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Figure 16. Similar to Fig. 1141 but the mass donors are TPAGB 

stars. 
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Figure 17. Mass ratio versus the mass for the BSs resulting from 
TPAGB stars transferring matter to MS companions for various 
/3 when the final mass of the primary is M%{ = M c + 0.2Mq. The 
small dots are at the end of RLOF. 



The influence of (3 on the BS formation from 
TPAGB+MS binaries is non-monotonic,i.e. the BS number 
first increases, and then decreases with the increase of j3 
(see Table [8]). This non-monotonic phenomenon is induced 
by the criterion of q c and the mass of the MS companion 
after accretion. With the increase of f3, q c becomes smaller 
and smaller, which leads q < q c to be fulfilled more un- 
likely, but the MS companion may accrete more material 
from the primary to increase the mass to be larger than the 
turn-off of the cluster. The results of Table |S] indicate that, 
when j3 < 0.4, the mass of the MS companion is the main 
factor to determine the formation of BS from this way, but 
the dynamical instability is the crucial factor when j3 > 0.4. 

From Fig. [I4]we see that, for the BSs from mass trans- 
fer between an FGB star and a MS companion, their masses 
increase and the range of the orbital period becomes nar- 
rower (mainly some long orbital period BSs disappear) with 
the increasing (3, since a large /3 means more matter ac- 
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creted by the secondary (leading to a larger final mass) and 
less angular momentum lost from the system (resulting in 
a shorter orbital period). Furthermore, a long-orbital pe- 
riod indicates a larger stellar radius at the onset of RLOF, 
and then a smaller q c as shown in section 3. Therefore, with 
the increase of /3, the RLOF in relatively long-orbital-period 
binaries firstly changes from dynamically stable to dynam- 
ically unstable, and directly leads to the disappearance of 
long-orbital-period BSs. The orbital period of the BSs from 
this channel has a wide range, i.e. from several days to hun- 
dreds of days for /3 — 0.1. However, BSs with P > 100 days 
can only be produced when /3 < 0.2 and their masses are 
near the turnoff because only a very small fraction of the 
lost matter from the primary is accreted by the secondary. 
All the possible candidates contributing to BSs in M67 from 
FGB+MS are low-mass binaries, since the secondaries from 
intermediate-mass FGB+MS binaries have left the main se- 
quence at 4 Gyr. For a 2Mq star, the minimum mass of the 
secondary is about 1.6M Q for (5 = 0.0 and M c = O.25M . 
The stars with this mass have left the main sequence in the 
old open cluster. 

Different from the case of FGB+MS, the BSs result- 
ing from EAGB+MS and TPAGB+MS have similar orbital 
periods for various /3 (see Figs. [15] and I16[l . Meanwhile, 
the mass difference of the BSs from different /3 is also not 
as obvious as that from FGB+MS binaries. In the cases 
of EAGB+MS and TPAGB+MS, some BSs have 
masses similar to those from FGB+MS, and some 
have lower masses, when /3 is larger than 0.1. The 
lower mass BSs mix with those from /3 with lower 
than the given value. Stellar wind may account for 
this consequence. When the primary evolves from 
FGB to EAGB, then to TPAGB, its mass becomes 
less and less due to stellar wind. Correspondingly, 
the minimum mass of the MS companion for stable 
mass transfer also becomes smaller. The companions 
may then have lower masses after RLOF with evo- 
lution. As a result, the masses of blue stragglers in 
Figs. 15 and 16 extend to lower values than those in 
Fig.14. 

From the criteria of the polytropic model and Hur- 
ley et al. (|2002ll . only long-orbital period BSs (with orbital 
period greater than 1600 d) may be produced from giants 
transferring matter to the MS companions. Since case A and 
early case B (the mass donor is during Hertzsprung gap at 
the onset of RLOF) mass transfer are only responsible for 
some short- and mid-period BSs (up to 100 d, Chen & Han, 
2004), there seems to be a period gap of hundreds days from 
the two criteria above. However, this period gap will not 
appear from the criterion of this paper. The BSs from differ- 
ent giant binaries cover the period range from several days 
to thousands of days (see Figs. H4l to [T6)l . 

As an example, we show the mass ratio versus the mass 
of BSs for the case of TPAGB+MS binaries at the onset and 
at the termination of RLOF in Fig. 1171 From the figure, we 
see that the mass ratio at the onset of RLOF is located in a 
triangle region, where the upper boundary is determined by 
the dynamical instability and the lower boundary is derter- 
mined by the mass of the accretor after the RLOF. 

We have not included example 3 in section 4 in the 
above study. Systematic investigation for this case is diffi- 
cult at present, since we have little knowledge of the condi- 



tions for binaries to avoid CE during initially dynamically 
unstable RLOF. The conditions are possibly relevant to the 
parameters such as the mass (and the core mass) of the pri- 
mary, the mass ratio of the components, the orbital period 
etc.. Furthermore, it may also be affected by the detailed 
process of mass transfer. The study of the conditions is be- 
yond the scope of this paper. If case 3 is a real case, only 
a very small fraction (less than 1 per cent) from those with 
dynamically unstable RLOF is enough to produce BSs in 
M67. The products here have much longer main-sequence 
lives in comparison to those from dynamically stable RLOF, 
since the secondaries in general have much lower mass due 
to the initially large mass ratio (greater than q c ) and are 
less evolved before RLOF. From Table [5j we see that at 
t = 2.47Gyr, a star with 1.5Mq has already formed and it 
becomes a BS when t > 2.5Gyr. l 13 IThe material around the 
BSs might give us some clues on their parent stars. 



5.4 Comparison to Observations 

Several works focus on the BSs in M67 from observa- 
tions (Milone & Latham 1992; Ahumada & Lapassct 1995 
|Sandquist fc Shetrone 2003 I. In the new catalogue of BSs 
in open clusters | |Ahumada fc Lap assct 2007}, there are 30 
BSs in this old open cluster. Since there is no orbital in- 
formation for these BSs in Ahumada & Lapasset i|2007[) . 
the main observational data we adopted here are from some 
earlier studies. 
Milone 



Latham (1992) reported the radial veloci- 
ties for 13 BSs in the open cluster M67 according to ob- 
servations over about 9 years. Three of the 13 BSs ro- 
tate too rapidly to allow reliable velocity determinations 
([Latham & Milone 1996). Among the ten BSs left, only one 
(F190) shows a short orbital period, about 4.2 d, while five 
have long-orbital periods in the range from 800 to 5000 d. 
The other four BSs are considered as single stars by some 
studies (e.g Hurley et al., 2001). Three of the five long-orbital 
period BSs have obvious orbital eccentricities and the other 
two are in near-circular orbits. Sandquist & Shetrone (2003) 
presented an analysis of the time series photometry of M67 
for W UMa systems, BSs and related objects. There are 24 
possible BSs observed in their study and most of them show 
no variation in their light curves. Two BSs, S968 and S1263, 
which were considered as single stars before, are possible 
variables in the study of Sandquist & Shetrone (2003). 

The BSs resulting from the criterion of Hurley et al. 
1(2002)) have orbital periods beyond 1600 d (logP > 3.2 as 
seen in Fig. I13[) . However, four of the five observed long- 
orbital period BSs have orbital periods less than this value. 
The simulation from the polytropic model gives a similar 
result. The BSs from the criterion of this paper may well 
cover the orbital period range, but it seems that none of the 
five BSs is from the mass transfer between an FGB and a 
MS, since the smallest orbital period is 846 d (logP = 2.93), 
which is far greater than the maximum orbital period result- 
ing from the FGB+MS binaries, as seen in Fig. [14] The BSs 
with mid orbital period (i.e. from tens of days to hundreds 
of days) are possibly from FGBs transferring matter to 

13 The age t = 2.5Gyr is for a cluster with a turnoff of around 
1.5M . 
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their MS companions, but it needs more observational 
evidence to comfirm the existence of these BSs. 

The orbital eccentricity is a puzzle from mass transfer, 
since the orbit should be circularized by tidal interaction of 
the two components of a binary, which is just undergoing, 
or after, RLOF. Some studies considered that dynamical 
collision is necessary to explain the observed orbital eccen- 
tricities. Recently, Marinovic et al. (|2007[> found that, due 
to the enhanced mass loss of the AGB component at orbital 
phases closer to the periastron, the net eccentricity growth 
rate in one orbit is comparable to the rate of tidal circular- 
isation in many cases. They reproduced the orbital period 
and eccentricity of the Sirius system with this eccentricity 
enhancing mechanism. Their study provides an explanation 
for the eccentricities of long-orbital period BSs without dy- 
namical collision. 



6 DISCUSSIONS AND CONCLUSIONS 

Our study shows that the critical mass ratio, q c , for dy- 
namically stable mass transfer between a giant star and a 
MS companion depends on the stellar radius at the onset 
of mass transfer. For any given mass, a more evolved star 
is less stable for RLOF, since the evolutionary timescale for 
a more evolved star is shorter and hence the mass transfer 
rate is higher than that for less-evolved ones. The re- 
sults including stellar wind ()Han et al. 2002[) are consistent 
with the tendency in this paper, except for cases where 
Mn = 0.8Mq, which have a larger / c and a larger q c in 
comparison to those of Mu > I.OMq. The behaviour of the 
cases of Mu = 0.8Mq here is similar to the consequences of 
a polytropic model. Meanwhile, from Table 3 in the paper 
of Han et al. (]2002| , we see that, at a similar core mass, bi- 
naries with a low primary mass (or a large f c )have larger 
q c , also consistent with equation(l) except for the binaries 
with Mu = 1.9Mq. Our results for j3 — in this paper 
match those of Han et al. (j2002f) (see the upper panel of 
Fig[5J. Due to the various core masses, however, the depen- 
dance of q c on f c is completely scattered as shown in the 
upper panel of Fig[9] So, it is very likely that both M c and 
/ c have influences on q c , but the effect is different for bina- 
ries of different masses, i.e. M c dominates the case where 
the primary's mass is larger than 1M© while f c is critical 
for the less massive ones. The thickness of the envelope 
might be an important cause here and the transition is pos- 
sibly a gradual process. The non-monotonic behaviour of 
Mu — 1.9M@ results from the low degeneracy degree of the 
core in primaries with this mass. 

From the criterion in this paper, we obtained 
some BSs with intermediate orbital period when j3 
is less than 0.4. However, there is no one reported 
at present located in the orbital period range from 
FGB+MS binaries, as we compare our results to ob- 
servations. There might be two factors for this con- 
tradiction. One is the fact that the mass of the BSs 
formed from such a low f3 is very close to the turnoff 
of M67, making them hard to distinguish from nor- 
mal stars around the turnoff (as seen in Figure 14). 
The other might be that the value of (3 is substan- 
tially larger than 0.1 during mass transfer between 
FGB stars and MS companions in a real case, so few 



BSs are produced in this way. Orbital determina- 
tions (in the future) might provide some constraints 
on the value of /3. For example, since the products 
of FGB+MS have a wide range of orbital periods, 
from several days to hundreds of days, there would 
be some BSs with intermediate orbital periods (if f3 
is low enough), which are possibly absent from the 
criterion of Hurley et al. (20021). 

The possibility that a main sequence star becomes a BS 
via wind accretion was first suggested by Williams (1964). 
This idea lacked attention for a long time, since an isolated 
star can hardly accrete enough matter to become a BS. How- 
ever it is likely different if the main sequence star is bound 
in a binary system where the primary is undergoing a large 
mass loss. In general, the mass loss rate of stellar wind is 
in the range of 10 _2 MQyr _1 to lO _6 M0yr _1 as a star ap- 
proaches the tip of the AGB, and some fraction of the lost 
material (up to 10 per cent, Tom et al. 1996) is probably 
accreted by its companion. This means that the secondary 
may significantly increase its mass when the primary under- 
goes a large mass loss at the tip of the AGB. Although the 
mass increase is likely not enough to produce a BS, the sec- 
ondary still has chances to obtain matter from the primary 
in the following stable RLOF if it happens. So BSs with very 
long orbital periods (greater than 1000 days) are likely the 
consequences of both RLOF and wind accretion. 

In this paper we present an example of a binary which 
avoids CE formation from initially dynamically unstable 
RLOF and eventually evolves to a BS with a long orbital 
period. We only assume that a CE has not formed in the 
initially dynamically unstable RLOF, but no critical con- 
ditions are shown. How to discriminate whether a CE has 
developed or not during this phase is unclear, as mentioned 
in section 5. Many parameters, such as the mass (and the 
core mass) of the primary, the mass ratio of the components, 
the orbital period etc. are relevant to this condition, which 
is also affected by the detailed process of mass transfer. The 
study of this will be complex and difficult, while inter- 
esting. As long as the possibility of this case exists, only a 
very small fraction of systems with dynamically unstable 
RLOF may provide an important contribution to BSs. 

Mass transfer efficiency, /3, is an important parameter. 
Both the criterion of dynamical instability, q c , and the fi- 
nal mass of the accretor, M a (which are the two critical 
characters to determine the formation of a BS) are 
relevant to this parameter. From our simulation for M67, the 
peak of the contribution from FGB+MS binaries is proba- 
bly about j3 = 0.1, with which the BS mass is near 1.3Mq. 
With the increase of j3, the orbital period decreases and the 
BS mass increases, as was explained in section 5. However it 
is unclear what the value of /3 is in binaries. In general, 
P will be very small for RLOF in binaries with a giant mass 
donor and a compact companion (Han et al. 20Q2J- So M a 
will not be very large, similar to the turnoff mass. 

The detailed evolution calculations in the paper show 
that q c decreases with the stellar radii of the primaries at 
the onset of RLOF, except for cases at or near the base of 
the giant branch, where the core is not very degenerate and 
the envelope is not yet fully convective. Non-conservative 
assumptions will strongly affect q c while stellar wind before 
mass transfer has little influence on it. To conveniently use 
the result we give a fit of q c as a function of the stellar 
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radius of the primary at the onset of RLOF, and of the mass 
transfer efficiency during RLOF. Theoretically, dynamically 
stable mass transfer occurs once the mass ratio is less than 
the critical value. However, it is delayed in real binaries. 
Usually the stable mass transfer occurs after the reversion 
of the mass ratio. 

The Monte Carlo simulations show that some binaries 
with the mass donor on the first giant branch, which have 
no contributions to the blue stragglers from the earlier cri- 
teria, will contribute to this population with the criterion 
obtained in this paper. Meanwhile, from our criterion, the 
blue stragglers resulting from the mass transfer between an 
AGB star and a MS companion may be more numerous and 
have a wider range of orbital periods than those formed 
from previous criteria. Although the result from our cri- 
terion may well cover the observed orbital period range, it 
seems that none of the five observed long-orbital period BSs 
is from the mass transfer between an FGB and a MS com- 
panion. 
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